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Abstract 

The interface between domains of opposite magnetization in the 3D Ising model 
near the critical temperature displays universal finite-size effects which can be de- 
scribed in terms of a gaussian model of capillary waves. It turns out that these finite- 
size corrections depend rather strongly on the shape of the lattice. This prediction, 
which has no adjustable parameters, is tested and accurately verified for various lattice 
shapes by means of numerical simulations with a cluster algorithm. This supports also 
a long-standing conjecture on the finite-size effects in Wilson loops of Lattice Gauge 
Theories. 
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1. Introduction 

The properties of interfaces in 3d statistical systems have been a long standing subject 
of research. In these last months there has been a renewed interest on the subject, with a 
particular attention to the case of the 3d Ising model, partially due to the improvement 
of Monte Carlo simulations: cluster algorithms and multicanonical simulations allow now 
to compare theoretical prediction on various physical quantities with reliable and precise 
numerical estimates [1-6]. In particular, the interface above the roughening temperature 
behaves like a free, fluctuating surface with a fixed topology determined by the boundary 
conditions, so these new algorithms provide us with a powerful tool for testing our ideas 
on the behaviour of fluctuating surfaces. 

Another reason of interest on this subject is due to the fact that in three dimensions 
ordinary spin models are related through duality to lattice gauge theories (LGT). In 
particular the physics of interfaces becomes in (LGT)'s the physics of Wilson loops, and 
the string tension of the 3d Ising gauge model @, |§ should have the same value of the 
interface tension of the 3d Ising spin model, as it will be confirmed below. 

In this letter we study the free-energy of the interface as a function of its shape. We 
will show that, using the capillary waves approach and the semiclassical analysis of ref. ||, 
it is possible to infer a functional form of the interface free energy for asymmetric lattices, 
which turns out to be in remarkable agreement with numerical simulations. We will also 
discuss implications of this result in the context of LGT's. 

This letter is organized as follows: in sect. 2 we discuss the capillary wave approach in 
the case of interfaces in asymmetric lattices. In sect. 3 we describe the simulation and the 
numerical results. In sect. 4 we use duality to make contact with the 3d Z2 gauge model 
and discuss the scaling behaviour of our results, while the last section is devoted to some 
concluding remarks. 
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2. Interface tension for asymmetric lattices 

The interface between two equilibrium phases above the roughening temperature is 
essentially a critical system: it may be freely translated through the medium and the long- 
wavelength, transverse fluctuations in the position interface, i.e. capillary waves, have a 
small cost in energy (hence cannot be neglected in calculations). They can be viewed as the 
Goldstone modes associated with the spontaneous breaking of the transverse translational 
invariance juj. 

The starting point of the capillary wave approach is the assumption that these fluctu- 
ations are described by an effective Hamiltonian proportional to the change they produce 
in area of the interface 



H/ksT = a I dx\ I dx 2 



k ctei/ \ox 2 



(1) 



'0 J 

where the field h(xi,X2) describes the interface displacement from the equilibrium position 
as a function of the longitudinal coordinates x\ and x 2 , L\ (L 2 ) is the size of the lattice 
in the x\ (x 2 ) direction and a is the (reduced) interface tension Q 

The interface free-energy increment due to interface fluctuations (which we shall term 
from now on capillary wave contribution to the interface tension) is then given by 

F cw = -k B Tlogtr e- H ' kBT , (2) 

The Hamiltonian of eq. (jl]) is too difficult to be handled exactly. However a crucial obser- 
vation is that this theory can be expanded in the adimensional parameter [oL\L 2 )~ x and 
the leading order term is the gaussian model, which turns out to be an accurate approx- 
imation on the range of parameters used in current lattice simulations. Then we replace 
eq.(|l|) with the oL\L 2 — > oo limit H — * Hq 



H G /k B T = - I dxi dx 2 
£ Jo Jo 



dh\ 2 / 01, ^ 
dxi) \dx 2 



(3) 



: Note that eq.(|l]) coincides with the Nambu-Goto string action in a special frame. 
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Within this approximation, the integration over h implied in eq. (|2|) can be done exactly. 
The integral is divergent but can be regularized using, for instance, a suitable generaliza- 
tion of the Riemann ^-function regularization (see e.g. ]ll|). The answer depends only 
on the geometrical properties of the boundary. The case we are interested in is that of 
a rectangle of lengths Li,L 2 with opposite edges identified (hence a torus) due to the 
periodic boundary conditions on the lattice. In this case the capillary wave contribution 
turns out to be 

F cw (L!,L 2 )/k B T = log (V^\ V (t)\ 2 ) +c ; T = ij- , (4) 

where r] denotes the Dedekind eta function: 

00 

T,(T)=q*l[(l-<n ; q = e 2 ™ , (5 ) 

n=l 

c is an undetermined constant, and we have assumed, without loss of generality, L\ > L 2 . 
Thus, the interface partition function Z(L\,L 2 ) may be written in the form 

Z(L h L 2 ) =e- aLlL *Z cw , (6) 
where the capillary wave contribution Z cw = e - F cw/k B T can ^ e eX p anc [ ec j as follows 

It 2 

Z cw (LuL 2 ) = CJ^-q-^ (l + q + 2q 2 + 3q 3 + 5</ 4 + 7q 5 + llq 6 + . . .) (7) 

where C is an undetermined constant and q = exp(—2irj^). This short expansion is 
largely sufficient for an accurate numerical evaluation of this contribution in the range 
of parameters we shall analyze. Eq.(|^) is scale invariant, indeed it only depends on the 
asymmetry parameter L\/L 2 and not on the absolute value of L\ or L 2 . As a consequence, 
it cannot be observed (since it only gives a constant contribution) if one measures the 
interface tension on square lattices of different sizes, as it was done in ref. (2], ||, ||, |6|, whose 
results are not affected by our calculation (and will be used in our following analysis). 

We want to stress that eq.(Q) (or at least its dominant contribution) is a rather well 
known result both in conformal field theories and in lattice gauge theories, where it was 



first discussed by Liischer, Symanzik and Weisz [|T^] as a finite-size correction term for 
Wilson loops (which, as we mentioned above, can be considered as the dual problem of 
interfaces, see e.g. pl| ). To the best of our knowledge, the first time it appeared in the 
context of interface physics is in [pHy , where the leading term of eq.(Q) was calculated 
using the Epstein's function regularization. For a comparison with numerical simulations 
we need the full form of eq. (Q) . 

Let us conclude this section summarizing the various assumptions and approximations 
involved in the capillary wave approach described above. First, we are assuming that 
the surface separating the two phases has a well defined localized shape. Second, we are 
neglecting bubbles and overhangs, namely we are assuming that h(x\,X2) is single- valued. 
Third, we are taking small fluctuations. Fourth, we are assuming that the h field does not 
self- interact. The latter two approximations can be controlled (at least in principle) by 
taking into account further orders in the expansion parameter (oL\L2)~ 1 . 

3. Monte Carlo simulation 

We studied the 3d Ising model on an ordinary cubic lattice of size L\ x L2 x L3 and 
periodic boundary conditions in all directions. We performed our Monte Carlo simulation 



using the Swendsen-Wang cluster algorithm [17]. L2 and L3 were kept fixed to the values 



L2 = 10 and L3 = 120 respectively, while L\ ranged from L\ = 18 to L\ = 32. As we shall 
see below these sizes are constrained by the requirement L2 < L\ <C £3, and the chosen 
values of L2 and L3 were the largest we could study with the CPU time we had at our 
disposal. We shall call in the following L3 the "time " direction. 
The Ising action is 

S = -£V(2)o-(£ + /i) (8) 
where the sum is over all sites x and all directions \i of the lattice, and the spin variables 
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can take the values cr(x) = dbl. The corresponding partition function is 

Z= E z- 0S ■ 0) 

<r(£)=±l 

The critical temperature of the model is estimated to be f3 c = 0.221654(6) fHf] . We studied 
the model in the broken symmetry phase, at (5 = 0.2275. This choice allowed us to make 
a direct comparison with the data of 0]. It is also well inside the scaling region of the 
dual gauge model, as discussed in 0, 

We evaluated the interface tension following the lines of ref . || (see also [l4j] ) looking 
at the energy splitting Eq due to the tunneling between the two vacua corresponding to 
the two opposite magnetizations of the system. This energy splitting can be measured 
looking at the time-slice correlations ( for this reason L3 has been chosen much larger 
than L\ and L2). 

More precisely we define 

L1L2 

where xi,X2,t are the coordinates in the L\,L2,L^ directions respectively. We extracted 
from our simulations the averages (SoSt)- Their t-dependence is given by 

(S S t ) = c [e- tEo + e -( L 3-^o] + Cl [ e -tEi + e -(i 3 -t)Ei] + . . . ( U) 

where the second exponent in each bracket is due to the periodic boundary conditions in 
the L3 direction, and we have explicitly written also the next to leading energy Ei, since 
it also turns out to be important in the range of data we measure. 

Finally the interface tension can be extracted from Eq by looking at its volume depen- 
dence In the case of a square section L\ = L2 = L this is given simply by: 

E (L) = Ce-° L2 . (12) 

This was the formula used in , where a remarkable agreement with numerical data was 
found. 



Si t^t: E E ^1^2, t) , (10) 

Xl = l 22=1 



Following the analysis of the previous section, the volume dependence of Eq on an 
asymmetric lattice (at the first order in the capillary wave approximation) is given by 



E Q (L U L 2 ) = Z cw (LuL 2 ) e- CTLlL2 (13) 

where Z cw was defined in eq.(|7|). Since Z cw is scale invariant, and simply reduces to 
a L-independent constant if L\ = L 2 = L, eq.(^) is completely consistent with eq.(|l2|), 
with the same value of the interface tension a. 

A simple way to check this picture is to evaluate the ratios R = Eq(Li, L 2 )/Eq(L, L), 
with L = \JL\li2- We get 

R{Ll,L2) = " VWlJlV) ) ■ (14) 



The values of Eq(L,L) are obtained extrapolating eq, (|l2|) , taking for C and a the values 
quoted in §, namely a = 0.0150(1) and C = 0.188(4). 

The ratios R, for various values of L\/L 2 , evaluated theoretically (Rxh) an d by Monte 
Carlo simulations (Rmc), are presented in tab.l. Error bars on Eq were estimated with an 
ordinary jacknife procedure. The same data are plotted in fig.l. The agreement between 
numerical data and theoretical predictions is rather impressive: note that the latter has 



no adjustable parameter, as it is shown by the last member of eq.(14). 

Let us mention few non-trivial features of the simulation. The approach outlined above 
allows very precise estimates of Eq and (with a lower precision) Ei, but, on the other 
side, time slice correlations are affected by very strong cross-correlations. These must 
be distinguished by the ordinary correlation in the Monte Carlo time, which are almost 
completely eliminated by the use of the cluster algorithm. Usually, cross-correlations are 
taken into account by weighting the data in the fit with the inverted cross-correlation 
matrix, but in this case cross-correlations were so strong that this method was almost 
useless. We decided to eliminate cross-correlations directly in the simulation by scattering 
the evaluation of time slice correlations in Monte Carlo time (a method we had already 



used in the context of the Ising gauge model ]?], H, where similar, even if less severe, 
problems occurred). 

Let us conclude this section with two remarks. 

First, it is possible to show that the constant Co in eq.([iT|) is related to the spontaneous 
magnetization. We checked that this relation is fulfilled in our samples. Moreover, fitting 
our time slices correlations with eq.(|Tl]), we always obtained the confidence levels between 
50% and 80%. Thus we are rather confident on the reliability of our estimates of Eq. 

Second, it is quite interesting to see how important is the contribution of the next to 
leading terms to the partition function Z cw . The difference between the contribution of 
the whole Z cw given in eq.(^) (full line in fig.l) and of its leading term q~ 12 (dashed line 
in fig.l) is much bigger than our error bars. Taking into account only this latter one would 
have a complete disagreement with the data and would erroneously reject the capillary 
wave picture. The greatest part of this difference is due to the square root term of eq.(|7]) 
(which is generated by the zero modes) . This lesson is particularly important in the dual 
context of Lattice Gauge Theories, where it is a common habit to take into account (if 
any) only the leading term of the Liischer-Symanzik-Weisz contribution. It could well be 
that some of the scaling violations observed in recent simulations of LGT could be due to 
this 0. 

4. Scaling and duality 

The 3d Ising model is related through duality to the 3d Z2 gauge model. Let us call 
(3 = -j^-j; the inverse temperature of the gauge model. (5 and (3 are related by 

The critical point (3 C = 0.2216 • • • is mapped by duality in C = 0.7614 

As we mentioned above, the capillary wave corrections have their parallel in the dual 

2 Indeed it has been shown in Q |^ that the string tension extracted by Wilson loops (Polyakov lines) 
taking into account the whole string corrections shows, as expected, a good scaling behaviour. 
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context of lattice gauge theory in the so called string corrections, namely the sum over all 
the possible fluctuations of the flux tube (thought as a string-like object) joining a quark- 
antiquark pair. Roughly speaking, the time world-lines of a quark-antiquark pair form a 
Wilson loop (or a pair Polyakov lines in the case of finite temperature LGT's); then, in 
order to evaluate the string corrections one has to sum over all the possible fluctuations of 
the surface bordered by the Wilson loop (or by the Polyakov lines). It is thus apparent that 
the only difference between the present analysis and those performed in LGT's is in the 
boundary conditions. In the case of Wilson loops one has fixed boundary conditions, hence 
a rectangular geometry. In the case of Polyakov line correlations one has fixed boundary 
conditions in one direction and periodic in the other, hence a cylindric geometry. Finally 
in the present analysis we have a toroidal geometry. In all these cases the asymptotic 
string (or capillary wave) contributions can be evaluated exactly (we refer to Q and to || 
for the discussion of the rectangular and cylindric cases respectively). 

In [[?J we evaluated the string tension of the Ising gauge model for various values of (3 
looking at the expectation value of Wilson loops. These data are collected in tab. 2 and 
plotted in fig. 2 together with those on the interface tension taken from ref.s [||, ||, ||. One 
can see that the two sets of data are perfectly compatible. The only difference is that 
those coming from the gauge model have bigger errors because the algorithms recently 
introduced for the calculation of the interface tension allow, with the same CPU time, 
an impressive improvement in precision. Notice that in the data taken from |7|, string 
corrections were already taken into account, hence the remarkable numerical coincidence 
of string and interface tension yields further support to the reliability of the capillary 
wave (or string in LGT's) approximation. Fitting the data of ref. Q according to the 
law cr(/3) = <r (l - (3 c /(3) 2u , with the fixed value of v = 0.63 |J|, we obtain a = 1.52(2) 
(x 2 /d.o.f. ~ 0.8) in perfect agreement with the results of [1-6]. 

It has to be noticed that in our analysis of Wilson loops |7L B ^(J we used a modified 
form of the gaussian model, where the field h is compactified on a circle of suitable radius. 
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As a consequence, the partition function has a form very different from that given in 
eq. ([?]). Surprisingly enough, in the range 1 < L\/L% < 5, the values of the ratios defined 
in eq.(|l4|) are almost the same of those calculated with eq.(|7]). 

5. Conclusions 

In this letter we have shown that the finite-size effects in the interface of 3d Ising model 
is well described by the capillary wave theory. We have shown that a simple, precise way 
to check this description is by looking at the interface tension on asymmetric lattices. The 
agreement between numerical data and theoretical prediction is quite remarkable once the 
whole capillary wave contribution has been taken into account. In particular we have 
shown that a relevant role is played by the subdominant logarithmic correction contained 
in F cw . 

It would be interesting to obtain our results with other simulation techniques like 
the multi-magnetic approach j| or the cluster update of boundary conditions proposed 
m §. In this way one would probably obtain higher precision and it would be possible to 
see higher order corrections. These could be calculated perturbatively by expanding the 
theory in the parameter (crLii^) -1 . The first correction is the quartic term coming from 
the expansion of the square root in eq.(||) (see e.g. 

It would be also interesting to repeat the same analysis with other models to see how 
universal are the predictions of the capillary wave approximation. Indeed it is well known 
that the roughening phenomenon, which can be described within the same approach, shows 
an high degree of universality. 
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Figure Captions 

Fig.l) The values of the ratios R defined by eq,(|i~4]) are plotted: the continuous line and 
the dotted line being respectively calculated using eq.(|7|) and taking into account 
only the g~T2 term of the same equation. The ratios evaluated by Monte Carlo 
simulations are also plotted with error bars. Note that the theoretical prediction has 
no free parameters. 

Fig. 2) The interface (string) tension is plotted for various values of /3. The data with 
error bars refer to the string tension of Wilson loops in the Ising gauge model and 
are taken from f7|. The other data refer to the interface tension in the Ising spin 
model and are taken from ref.s 0, |3|, ^| , their error bars are not reported since are 
smaller than the plotted symbols. The continuous line is the fit of the data of ref |?j 
only with the law a(/3) = a Q (l - &//3) 1 - 26 . 

Table Captions 

Tab. I The same data of fig.l. are reported: in the first column the ratio L1/L2, in the 
third and fourth columns the corresponding ratios from the Monte Carlo (Rmc) an d 
theoretical (RTh) evaluation, while the values obtained taking into account only the 
q~T2 term is given in the last column. In the second column are written the values of 
the energy splitting Eq that we measured on the corresponding asymmetric lattices. 

Tab. II The same data of fig. 2. are reported: in the first column the values of (3, in the 
second the corresponding values of the interface (string) tension and in the third the 
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values of do obtained using the relation uq = cr(/3)(l — /? C //?)~ L26 . In the last column 
are reported the ref.s from which the data are taken and the simulation techniques 
adopted: WL=Wilson loops, MM=Multimagnetic simulation, CUBC=Cluster up- 
date of boundary conditions, ES=energy splitting. 



Tab. I 



x = L 1 /L 2 


E (L U L 2 ) 


Rmc(x) 


R T h(x) 




18/10 


0.01439 (26) 


1.14 (6) 


1.13 


1.52 


20/10 


0.01099 (23) 


1.17 (7) 


1.19 


1.69 


23/10 


0.00797 (23) 


1.33 (10) 


1.30 


1.97 


26/10 


0.00584 (23) 


1.53 (13) 


1.43 


2.31 


28/10 


0.00396 (44) 


1.40 (22) 


1.53 


2.57 


32/10 


0.00289 (32) 


1.87 (31) 


1.76 


3.16 



Tab. II 



(3 


a 


O"0 


meth. & ref. 


0.2230 


0.00252 (3) 


1.586 (19) 


CUBC [6] 


0.2240 


0.00492 (4) 


1.529 (12) 


CUBC [6] 


0.2255 


0.00904 (6) 


1.529 (10) 


CUBC [6] 


0.2258 


0.0093 (5) 


1.442 (78) 


WL [7] 


0.2269 


0.0135 (5) 


1.538 (57) 


WL [7] 


0.2270 


0.01293 (17) 


1.457 (19) 


MM [3] 


0.2275 


0.0150 (1) 


1.508 (10) 


ES [2] 


0.2281 


0.0177 (5) 


1.575 (45) 


WL [7] 


0.2293 


0.0212 (5) 


1.537 (37) 


WL [7] 


0.2305 


0.0244 (5) 


1.484 (31) 


WL [7] 


0.2317 


0.0292 (6) 


1.523 (37) 


WL [7] 


0.2320 


0.03140 (14) 


1.575 (7) 


MM [3] 


0.2327 


0.0328 (2) 


1.525 (10) 


ES [2] 
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